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Abstract 

Using the pseudospin representation and the SU(2) phase operators we introduce a complex 
parameter to characterize both infinite and finite superconducting systems. While in the bulk 
limit the parameter becomes identical to the conventional order parameter, in the nanoscopic limit 
its modulus reduces to the number parity effect parameter and its phase takes discrete values. We 
evaluate the Josephson coupling energy and show that in bulk superconductor it reproduces the 
conventional expression and in the nanoscopic limit it leads to quantized Josephson effect. Finally, 
we study the phase flow or dual Josephson effect in a superconductor with fixed number of electrons. 

PACS numbers: 74.20.Fg, 74.20.-z, 74.50. +r 



Recent experimental works on superconducting metallic islands at nanometer scale have 
established a link between bulk superconductors and atomic nuclei as regards pairing 
correlations^**- The long-range order in a bulk superconductor can be described by an 
order parameter A which is complex and the equations have symmetry properties which 
ensure that if A is a solution, then e t9 A is also a solution^. On the other hand when size 
of a superconductor is reduced to nanometer scales so that the number of the electrons is 
fixed, the order parameter A = (c_k|Ck|) vanishes where c_kj. and Ckf are the annihilation 
operators for time-reversed states | — k |) and | k j), respectively. In this case, supercon- 
ductivity manifests itself with nonvanishing number parity effect parameter Ap where the 
ground state energy of the system increases or decreases, depending upon whether the total 
number becomes odd or even, by addition of a new electron^. In this work, we propose 
a parameter, which unifies the order parameter A of the bulk limit and the number par- 
ity effect parameter Ap of the nanoscopic superconductors. Introducing the pseudospin or 
quasi-spin representation for the model Hamiltonian of the theory of superconductivity^ 
and the SU(2) phase states^ we define a quantum phase for a superconductor with discrete 
energy levels along with modulus of the parameter which becomes equal to Ap. As we 
go from the nanoscopic limit to the bulk superconductor we show that the number parity 
effect parameter and the SU(2) phase go to the amplitude and the phase of the bulk order 
parameter, respectively. 

We are going to start with a notation which is more proper for nanoscopic superconductors 
where energy levels are discrete and finite^. This reduced form of the BCS model was applied 
in nuclear physics and it has an exact solution 12 . The model Hamiltonian is 

J> 3>j' 

where g is the pairing coupling constant for the time-reversed states | j f) and | j j), both 
having the energy €j. Here, cJ CT (cj a ) is the creation (annihilation) operator for state | ja) 
where j G {1, Vl} and a E {T, !}• 

Introducing the pseudospin variables^ 

s j = \ (sW + c )i c n ~ l ) > s 7 = c n c tf = (4) t ( 2 ) 

which generate the SU(2) algebra 

= ±5ijsf, (3) 
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it is possible to rewrite the model Hamiltonian as 

H = ^U+ 1 -)-gY,4sJ- (4) 

We note that the mapping from the Fermi operators to the pseudospin operators is possible as 
long as all single particle states are doubly occupied. However, since the original Hamiltonian 
(JTJ contains no terms which couple a singly occupied level to others, the only role of such 
states will be blocking from pairing interaction. Therefore, the summations in Eqn. (J3| are 
over doubly occupied or empty states. Both the above mapping and the BCS wave function 4 
lack proper antisymmetrization due to separate treatment of singly occupied states, but 
since the model Hamiltonian (0) does not involve any scatterings into or out of such states, 
antisymmetrization with respect to interlevel pair exchange and intrapair electron exchange 
is sufficient. 

In this work, rather than the exact solution of the problem, we are interested in the 
qualitative result which has also been obtained numerically: The ground state energy for 
even number of electrons is lower in comparison to neighboring odd number state o 13 i 14 i 15 i 16 
including degenerate case^-. Parity dependence of the condensation energy and pairing 
parameters in nanoscopic superconductors was first emphasized by von Delft et a/.— but the 
first correction to the bulk limit had been obtained by Janko, Smith and Ambegaokar— and 
Golubev and Zaikin 19 . 

Phase operators and phase states have been studied mainly in quantum optics and possi- 
ble connection of quantum phase and the mean field treatment of the BCS Hamiltonian has 
been pointed out by Shumovsky 2 -. Given SU(2) algebra, for example the one generated by 
the components of the total spin operator s = J2i s i, we can introduce^ the radial operator 
defined by 

s r = V s + s~ (5) 
and the exponential of the phase operator given by 

m=s 

E= I S;sm + l)(S;sm | . (6) 

m=—s 

Here, | S; sm) is simultaneous eigenstate of s 2 and s z operators with eigenvalues s(s + 1) 
and m, respectively. In order to simplify the notation, m is defined modulo 2s + 1 so that 
| S; ss + 1) =| S; s — s). The label S has been introduced to distinguish them from the phase 
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states to be defined below. We are going to make use of the cosine and the sine operators 

Cos =\{E + E- 1 ) , Sm=^(E- E- 1 ) . (7) 

For integer s or on the so called Bose sector, the eigenstate of E with eigenvalue 
exp(— i27r/x/(2s + 1)) is evaluated to be 



I m=s 

9; su) = ■ exp 



2tt/j 

i m 

2s + 1 



S; sm) (8) 



V2s~T 

and a similar expression holds for half integer s or in the Fermi sector. 

In terms of the radial and the exponential of the phase operators for the total spin, it 
is possible to rewrite the interaction part of the Hamiltonian (0} as —gs r EE^s r . Since E 
is unitary, we obviously have EE^ = I but our aim in keeping E and E^ is to define the 
phase properly. Now, we introduce the superconductivity criterion as (s r ) ^ 0. We are 
going to prove that this definition agrees with existing criteria for both grand canonical and 
canonical superconducting systems. We are going to show that (s r ) becomes identical to 
the modulus of the BCS order parameter in the bulk limit while in the nanoscopic limit it 
reduces to the number parity effect parameter Ap in units of g. There have been several 
suggestions for a canonically meaningful pairing paramete r 11 ' 1 ^ 21 ' 22 ' 2 ?' 24 . Our definition is 
equivalent to that of Tian et. al^, which has been proposed by Penrose and Onsagei^ and 
Yang2& as a measure of the strength of the spontaneous symmetry breaking field. Amico 
and Osterloh 2 ^. and Zhou et.a® 8 have calculated the pairing correlation function (s7"sj) 
analytically by extending Richardon's results^. We further introduce (E) as exponential 
of the phase. This definition is justified by the observation that in the grand canonical 
ensemble (E) turns out to be exponential of phase of the BCS order parameter. 

We first note that [s r , s z ] = and hence s r gives a good quantum number even for a 
finite system. Secondly, (s r ) is filling dependent even for a single, d— fold degenerate level in 



contrast to Ap. While (s r ) = Jv(d — v + 1) for v pairs, the number parity effect parameter 
Ap is gd/2, independent of v. For v — d/2, i.e. half-filling or m = 0, the two results become 
identical. 

We start to our proof by examining the canonical system. In analogy to the pairing 
energy in nuclear physics*^, Matveev and Larkir>£> introduced the parity effect parameter 

Ap = E 2n+l - E2n \ E2n+2 (9) 



for nanoscopic superconductors where E n is the ground state energy for n electrons. Assum- 
ing that the expectation value of the single particle energy part (the first term in Eqn. (|3|)) 
follows a monotonic behavior so that T 2n+ i — (T 2n + T 2n+2 )/2, the main contribution to the 
ground state energy will come from the interaction part so that 

Ap = -g I (s r ) 2 n+i - I • (10) 

Now, the eigenstates and in particular the ground state of the model Hamiltonian will be of 
the form 

X)c sm | S] sm) (11) 

because the interaction term commutes with s 2 and s z while the single particle part com- 
mutes with the latter only and hence m is a good quantum number. We note that since s is 
the total spin in general it is multiply degenerate. It is possible to calculate the expectation 
value of the radial operator s r = \J s + s~ as 

(s r )n = J2\c sm | 2 <Js(s + 1) -m(m- 1). (12) 

s 

Here, the number of electrons n is a function of m. In BCS theory the single particle states 
participating in pairing interaction are assumed to be those in a shell of thickness ~ 2hu>D, 
uj£, being the Debye frequency, symmetric around the Fermi level. In this case, half of the 
states are full while the half is empty and hence m = 0. Near half-filling where m ~ and 
for s 1, we can approximate the square root as s to give (s r ) n ~ J2 S I c so | 2 s - Similarly, 
for s 2 , with the same approximations we find that 

(S 2 r )2n ^ (s 2 r }2n+2 ~ ^ I C *0 ^ ^ ■ ( 13 ) 

s 

For 2n + 1 electrons, the mere effect of the unpaired electron is to block one of the single- 
particle energy levels from pairing which in our notation means that the corresponding spin 
value becomes s — 1/2. However, using Eqn. ()10j) this simply gives that Ap ~ g(s r )- 
Therefore, the parameter we introduced (s r ) (multiplied by the pairing coupling constant g) 
is identical to the number parity effect in the proper limit. 

Next, we examine the grand canonical case or the thermodynamic limit. In its present 
form the model Hamiltonian (j3J) commutes with s z and therefore m is a good quantum 
number or equivalently the number of electrons is a conserved quantity. To make a con- 
nection with the BCS order parameter we are going to replace the interaction part of the 
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Hamiltonian by 



g(s r EA* + AEls r - | A | 2 ) 



(14) 



which is nothing but the standard mean field approximation since s r E = s + and E + s r = s~. 
The BCS wave function describes a state with totally indefinite number of particles but with 
a definite phase. We can project the BCS states onto states of definite particle number by 
taking the Fourier transform with respect to the phased and that is why particle number 
iV and phase are conjugate variables with an uncertainty relation SNS(p ~ 1. It has been 
shown that in the thermodynamic limit the ground state of the BCS Hamiltonian (j^J) is 
also the ground state of the mean field Hamiltonian whose interaction part is given by Eqn. 
(J14|) 32 . These are nothing but the phase states which we have defined above. In our case this 
result can be verified by observing that near half-filling and at large s, we have [s r , E] ~ 0. 
Therefore, we evaluate the expectation value of s r E in state | 9; s/j) and find that 



We identify the phase —2ii[i/ (2s + 1) as <fi and the factor in front (the sum divided by 2s + 1) 
as the modulus of the order parameter | A | . This completes our argument on the relation 
of (s r ) and (E) to Ap and A except one point: What happens to {E) for a system with 
discrete energy levels but yet with indefinite number of electrons? We note that this not 
the thermodynamic limit. The system is finite but yet the number of electrons is not fixed. 
Such a situation can be realized through a Josephson junction. 

The origin of the Josephson interaction is single-particle tunneling electron pairs. At low 
energies, single-particle tunneling interaction lead to two contributions both of which are 
second order processes. The first one, where an electron goes from one superconductor to 
the other and returns, leads to proximity effect. The second one is the Josephson tunneling 
of two electrons from one superconductor to the other. The only effect of the first process is 
to renormalize the single-particle energies. Furthermore, there is no net current associated 
with this process. We can evaluate the explicit contributions of these two processes by 
considering two superconductors, both of which are described by the model Hamiltonian Q 
so that we are going to denote the total Hamiltonian as Hq. Let us consider a tunneling 
interaction of the form 




(15) 



E ( C 2j>a C lj° + C lja°2j'a) 



(16) 



where Cij a (c2j a ) is the annihilation operator for state | ja) of the first (second) superconduc- 
tor. One way to introduce V perturbatively is to use the unitary transformation method- 3 
where the second order Hamiltonian takes the form H + [V, f2]/2. Here, the anti-Hermitian 
operator Q is given by 

^ | mim 2 )(mim 2 | V | mn 2 )(nira 2 | , , 

^~ 2^ (o) (o) 

where | nin 2 ) denotes the ground state of H with ri\ electrons in the first superconductor 
and n 2 electrons in the second and e^ n is the corresponding energy eigenvalue of the 
combined system. Since we are interested in low energy excitations, at each step we project 
the system into its ground state. The two contributions we discussed above, the proximity 
and Josephson processes, can easily be calculated. Repeating the approximations we did in 
Eqn. (JT3J), we find that the strength of both terms are given by — t 2 /(Aip + A 2 p) = ej. 
In particular the Josephson interaction term can then be written as Sji^EiE^ 1 + Ei 1 E 2 ))/2 
where Ei is the exponential of the phase operator in the i th superconductor. We immediately 
observe that for phase state | 0i0 2 ), expectation value of this term is simply e j cos(0i — 2 ). 
To simplify our final analysis let us assume that one of the superconductors is large so that 
it can be described by the BCS state with a fixed phase which we can assume to be zero 
without loss of generality. Then for the other we can write down an effective Hamiltonian 

H eff = £ 2i J ( s j + \)~ 9s + s- + 4e c (s z - (s z )) 2 + SjCos (18) 

where Eq is single-electron charging energy of the island and ij is renormalized single-particle 
energy. The Josephson current Ij = 2e(s z ) can be easily calculated as Ij = 2esj(Sin) /Ti 
where Sin is the sine operator. Therefore, in the bulk limit where eigenstates are nearly 
phase states, we recover the conventional expression for the Josephson curren t? 4 ^ . Eigen- 
states of the Hamiltonian composed of the first three terms of H e ff are still given by (jll|) 
where m is a good quantum number. Hence, H e ff is nothing but tight-binding Hamiltonian 
with nearest neighbor (mil) hopping matrix element Sj/2. The nature of the eigenstates 
depends upon the on-site energies. For example, for quadratic dependence of energy eigen- 
values (in the absence of SjCos term) on m, which would be the case for flat ij, we can 
find the exact eigenvalues and eigenstates of H e ff. In this case we obtain a tight-binding 
Hamiltonian with on-site energies having quadratic dependence on site index and we can 
find the solution by observing that the expansion coefficients of the Mathieu function ce 2n 
satisfy a recursion relation which is identical to the characteristic equation of H e ff&. 



It is clear that Ij vanishes for any state which can be written as a linear combination 
of | S; sm) states with real expansion coefficients. These are nothing but bound states in 
S— space. On the other hand for propagating states, like | 9; sfi), Ij is non-zero. For small 
enough systems a very interesting situation may arise because discreteness of fi and hence 
quantized Ij values might be observed. In other words, if the number of the single particle 
energy levels and hence s is not too big, we can measure a quantized Josephson current. 
A single electron transistor with a small enough superconducting island can be used to see 
quantization effect. Another possibility is to measure the Josephson plasma oscillations 
between a bulk and nanoscopic superconducto r 37138 . Recent first principle calculations for 
structural and electronic properties of aluminum covered single wall carbon nanotubes show 
that a stable metallic ring can be formed 3 ^. These structures can also allow us to observe 
effect of phase quantization. Coulomb interaction works in the direction to suppress the 
current but using an external electric field relative strength of the Josephson interaction can 
be increased. One possibility is to measure the Josephson current through one dimensional 
array of aluminum rings formed around a carbon nanotube. In general, any physical quantity 
depending upon the phase is a candidate to observe quantization. For BCS gap A = 
2hurj£~ 1 ^ x and level spacing 5 satisfying 5 ^ A, assuming for example that we can resolve 
discreteness of the phase angle for 2Hujd/S ~ 1000 states in the Debye shell at the Fermi 
level, we evaluate A to be > 0.14. For larger A, we can go to smaller sizes or less number of 
states and hence there is more chance to observe quantization effects. 

If N and are conjugate variables and the Josephson effect is a phenomenon relevant to 
fixed (f) and indefinite N, what is its dual effect where N is fixed but <p is indefinite? When 
a nanoscopic superconductor is coupled by Coulomb interaction to another superconductor, 
there appears a second order effect which is analogue or dual of Josephson effect where 
particle numbers are fixed but the phases are not determined. To make the analogy complete 
let us consider an interaction term of the form Ed{F + F" 1 ) where Ed is the dual Josephson 
interaction energy and F is dual to the operator E and it is defined by 



F = exp i 



2s + 1 



(19) 



It is easy to show that F \ 9; sfi) =\ 9; sfi + l) 10 . We can evaluate the phase current (9 Z ) 



where 



s 



®z= ^ I 0; s/j)(9;s/i 



(20) 



fl= — S 
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in complete analogy to Josephson current as ieoiF^ 1 — F)/H. This interaction is similar to 
the van der Waals force between two molecules which is a manifestation of discreteness of 
electronic energy levels. In the superconducting state, intragrain single particle excitation 
spectra are modified due to the number parity effect and hence there appears an additional 
interaction due to pairing. In other words, dual Josephson effect refers the attractive in- 
teraction between two superconductors due to virtual Cooper pair breaking (as a result of 
Coulomb interaction between the superconductors) where interaction energy is of the 
order of the ratio of Coulomb interaction squared to superconducting gap or number parity 
effect parameter. This effect might also have relevance to atomic nuclei when they approach 
close enough so that Coulomb force becomes appreciable. 

In conclusion, we proposed a complex parameter to describe pairing correlations in a 
fermionic system. We showed that our definition agrees with the existing parameters in the 
canonical and grand canonical descriptions. We predicted possible quantization in Josephson 
effect in the nanoscopic limit. We further analyzed the dual Josephson effect a nanoscopic 
superconductor and interpreted the resulting expression in terms of quantum phase flow. 
Recently, the complex parameter introduced this work has been used to study quantum 
entanglement a paired finite Fermi system 4 ^. 
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